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Abstract 

We discuss kinematical correlations between charged leptons from semileptonic decays of open 
charm/bottom, leptons produced in the Drell-Yan mechanism as well as some other mechanisms 
not included so far in the literature in proton-proton scattering at BNL RHIC. The distributions 
of charm and bottom quarks / antiquarks are calculated in the framework of the /^-factorization ap- 
proach. For this calculation we use different unintegrated parton distributions from the literature. 
The hadronization of heavy quarks is done with the help of well-known fragmentation functions. 
Uncertainties of our predictions related to heavy quark masses, factorization and renormalization 
scales as well as due to the choice of fragmentation model are also discussed. We use semileptonic 
decay functions found by fitting recent semileptonic data obtained by the CLEO and BABAR 
collaborations. The Drell-Yan processes were calculated including transverse momenta of quarks 
and antiquarks, using the Kwiecihski parton distributions. We have also took into consideration 
reactions initiated by purely QED 7*7*-fusion in elastic and inelastic pp collisions as well as re- 
cently proposed diffractive mechanism of exclusive charm-anticharm production. The contribution 
of the later mechanism is rather small. We get good description of the dilepton invariant mass 
spectrum measured recently by the PHENIX collaboration and present predictions for the dilepton 
pair transverse momentum distribution as well as distribution in azimuthal angle between electron 
and positron. 
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I. INTRODUCTION 



Recently the PHENIX collaboration has measured dilepton invariant mass spectrum from 
to 8 GeV in proton-proton collisions at y/s = 200 GeV [1] . It is commonly believed that the 
main contribution to the dielectron continuum comes from so-called nonphotonic electrons 
which are produced mainly in semileptonic decays of charm and bottom mesons. Up to 
now, production of open charm and bottom was studied only in inclusive measurements of 
charmed mesons [2] and electrons [3] and only inclusive observables were calculated in pQCD 
approach [4, 5]. Such predictions give rather good description of the experimental data, 
however, the theoretical uncertainties are quite large which makes the situation somewhat 
clouded and prevents definite conclusions. 

Some time ago we have studied kinematical correlations of cc quarks [6], which is, however, 
difficult to study experimentally. High luminosity and in a consequence better statistics at 
present colliders gives a new possibility to study not only inclusive distributions but also 
correlations between outgoing particles (meson- meson, meson-electron or electron-electron). 
Kinematical correlations constitute an alternative method to pin down the cross section 
for charm and bottom production. It gives also a great possibility to separate charm and 
bottom contributions which has a crucial meaning for understanding the character of heavy 
quarks interactions with the matter created in high energy nuclear collisions [7]. 



II. FORMALISM 

The /^-factorization method is very useful to study correlations between cc [6] and e + e~ 
from the Drell-Yan processes [8]. In our calculations we take under consideration not only 
leptons from open charm/bottom decays but also leptons produced in Drell-Yan proccess, 
as well as leptons coming from elastic and inelastic processes initiated by photon-photon 
fusion. In the case of elastic reaction we follow exact momentum space calculations with 4- 
body phase space (see e.g. [9]) and for inelastic scattering we have applied unique (collinear) 
photon distributions in the nucleon MRST2004 [10]. 



A. Dileptons from semileptonic decays 

The electrons from semileptonic decays are produced in a three-stage process. The whole 
procedure can be written in the following schematic way: 

A?^<E>/^e, (2.1) 



dycPp dyd 2 p 

where the symbol <8> denotes a generic convolution. The first term is responsible for produc- 
tion of heavy quarks/ ant iquarks (see Fig.l). Next step is the process of formation of heavy 
mesons and the last ingredient is semileptonic decay of heavy mesons to electrons/positrons. 
The inclusive production of heavy quark/antiquark pairs can be calculated in the framework 
of the /^-factorization [11]. In this approach transverse momenta of initial partons are in- 
cluded and emission of gluons is encoded in a so-called unintegrated gluon, in general parton, 
distributions. In the leading-order approximation within the ^-factorization approach the 
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differential cross section for the QQ or Drell-Yan process can be written as: 

da sr^ f d 2 Ki <t d 2 K 2 ,t 1 



dy 1 dp lt dy 2 dp 
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where Ti{x\, K 2 t ) and J r j(x 2 , ft 2 t) are the unintegrated gluon (parton) distribution functions 
(UPDFs). The longitudinal momentum fractions can be calculated as 



x 2 



mu , , , m 2t , , 
-7=exp(yi) + -^exp(y 2 ) , 

m u m 2t , \ 

exp(-2/i) + — exp(-2/ 2 ) 
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where ?/i and y 2 are rapidities of heavy quark and heavy antiquark, and m u and m 2t are 
their transverse masses. 




FIG. 1: The dominant mechanism of the cc and bb production at high energy. The emission 
of several extra gluons is included in the unintegrated gluon (parton) distributions used in the 
formalism. 



There are two types of the LO 2 — > 2 subprocesses which contribute to heavy quarks 
production, gg — > QQ and qq — > QQ. The first mechanism dominates at large energies and 
the second one near the threshold. Only gg — > QQ mechanism is included here. We use 
off-shell matrix elements corresponding to off-shell kinematics so hard amplitude depends 
on transverse momenta (virtualities of initial gluons) in the exact way. At relatively low 
RHIC energies rather intermediate x- values become relevant so the Kwiecinski UGDFs seem 
applicable in this case [12]. However, to show the uncertainty of our predictions resulting 
from different approaches in calculating uninegrated parton distributions we have also used 
Kimber-Martin-Ryskin (KMR) [13] and Kutak-Stasto models [14]. All of them have different 
theoretical background. It is therefore very interesting to compare such results with the 
PHENIX data and verify applicability of these UGDFs at RHIC. In the case of the Kwiecinski 
distributions we fix the renormalization and factorization scales to standard values y? R = 
fi 2 F = Aitlq. Using the KMR UGDF the mostly used set of these parameters in the context of 
inclusive heavy quark production is fi 2 R = ^rrig and n 2 F = M^q, where Mqq is the invariant 
mass of the QQ pair. 
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The hadronization of heavy quarks is usually done with the help of fragmentation func- 
tions. The inclusive distributions of hadrons can be obtained through a convolution of 
inclusive distributions of heavy quarks/ant iquarks and Q — > h fragmentation functions: 



dcr(yi,Pi t ,y 2 ,P2t, ( P) 
dyidp^ t dy 2 dpgd<p 



I 



Dq^h{zi) _ Dq^h(z2) _ d<r(y 1 ,p < ? t ,y 2 ,P2t,<f>) 
z\ z 2 dy 1 dp^ t dy 2 dp2 t d(f) 



dz\dz\ 



■2 ■ 



(2.4) 



H 



where: pf t = p\\ = , where meson longitudinal fractions z±,Z2 G (0, 1). We have made 
approximation assuming that y±,y2,(f) are unchanged in the fragmentation process. 

There are several models of fragmentation functions in the literature. Here we mostly use 
the Peterson fragmentation function [15]. However, to check the sensitivity of our results 
to the choice of the fragmentation model, we have also applied fragmentation functions 
proposed by Kartvelishvili et al. [16] and Braaten et al. [17]. 

Recently the CLEO and BABAR collaborations have measured very precisely the spec- 
trum of electrons/positrons coming from the weak decays of D and B mesons, respectively 
[18]. These functions can in principle be calculated. This introduces, however, some model 
uncertainties and requires inclusion of all final state channels explicitly. An alternative is 
to use proper experimental input which after renormalizing to experimental branching frac- 
tions can be use to generate electrons/positrons in a Monte Carlo approach. The electrons 
(positrons) are generated isotropically in the heavy meson rest frame. In the present paper 
we use parametrizations of the decay functions found in Ref.[5]. 

B. Drell-Yan dileptons 

The electron and positron produced in the Drell-Yan mechanism are naturally correlated. 
We have shown recently [8] how to use the transverse momentum dependent parton (quark, 
antiquark) distributions to obtain several differential distributions. A basic diagram of 
the mechanism is shown in Fig. 2. In our calculations here we follow Ref.[8] and use the 
Kwiecinski parton distributions. Here the off-shellness of quark/antiquark is included in the 
kinematics and the matrix element taken here in the on-shell form expressed in terms of 
the subprocess invariants calculated with the off-shell condition. This is not fully consistent 
but avoids problems when on-shell momenta of quarks and on-shell matrix element are used 
[19, 20]. In any case the result of our approach is not very different than that for the 
collinear approach but includes kinematical effect of transverse momenta which is crucial 
to understand e.g. azimuthal correlations and distributions in transverse momentum of the 
dilepton pair, impossible to address in the collinear approach. 

The differential cross section for the 0-th order contribution including quark/antiquark 
transverse momenta can be written as: 



where T^Xi,^) and J r i{x 2 ,K\ t ) are unintegrated quark/antiquark distributions in hadron 
hi and /12, respectively. 
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The longitudinal momentum fractions are evaluated in terms of final lepton rapidities 
and transverse momenta: 

mu , v , m 2t , \ 
x 1 = -=exp(j/i) + —=exp{y 2 ) , 

m u , . m 2t , , 
£2 = - 7 =exp(-y 1 ) + — =exp(-y a J, 
V s V s 

where m t = a/j^ 2 + m 2 are transverse masses of electron and positron. 
The delta function in Eq.(2.5) can be eliminated as e.g. in Refs.[6]. 



(2.6) 




FIG. 2: The Drell-Yan mechanism of the dielectron pair production. The extra gluon emissions 
are included in the formalism of unintegrated parton distributions. 



C. QED 7*7* elastic and inelastic production of dileptons 

The matrix element for the pp — > ppe + e~ reaction via 7*7*-fusion (see Fig. 3) can be 
approximately written as 

« efiftO^^^foft) {P2+ / 2T eF ^ > ( 2 - 7 ) 

where F\{t\) and Fi{t 2 ) are Dirac proton electromagnetic form factors, and the 7*7* — > e + e~ 
vertex has the form 



ta7*7* 2 - n \ I v Qi - h - m qi - k 2 + m , « 



(2.8) 



The above formulae are used to calculate differential cross section via exact integration in 
the full 4-body phase-space. The details can be found e.g. in Ref. [9]. We shall call this 
contribution double elastic for brevity. 

In addition, there are components when one of the protons, or even both (see Fig. 3), 
do not survive the collision. Corresponding contributions will be called single and double 
inelastic processes respectively. The double inelastic contribution can be calculated as the 
gluon-gluon contribution in the parton model by replacing gluon distributions in the nucleon 
by corresponding photon distributions. Only one group discussed photon distributions in 
the nucleon [10]. The corresponding cross section can be calculated as 

da 



dyidy 2 d 2 pt 16tt 2 s ; 



;Xif 7 / P (xi,fi )x 2 / 7 / p (x 2 ,/i )\M^ e + e -\ 2 . (2.9) 
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FIG. 3: Diagrammatic representation of processes initiate by photon-photon subprocesses: double- 
elastic, double- inelastic, inelastic-elastic and elastic-inelastic. 



The matrix element can be found in several text books. The cross section for single 
inelastic process can be calculated by a replacement of one of / 7 / p (x,/i 2 ) by f^ l / p (x) which 
is often called elastic photon flux factor. Relevant formulae can be found in [21]. 

D. Exclusive double-diffractive production of open charm 

Recently, our group has calculated, for the first time in the literature, the exclusive double 
diffractive (EDD) production of open charm [22]. A sketch of this mechanism is shown in 



FIG. 4: The mechanism of exclusive double-diffractive production of open charm. 

The pp — > p(cc)p reaction is treated as a genuine 4-body process with exact kinematics. 
This can be easily used to apply kinematical cuts required by experiments. According 
to the Kaidalov-Khoze-Martin-Ryskin (KKMR) approach used previously for the exclusive 
Higgs boson production [23], the amplitude of the exclusive diffractive qq pair production 
pp — > p{qq)p can be written as [22] 



diagonal unintegrated gluon distributions in nucleon 1 and 2, respectively. 

The vertex factor V^ 2 = V^ 2 (qi,q2,ki,k 2 ) in expression (2.10) is the production am- 
plitude of a pair of massive quark q and antiquark q with helicities X q , Xq and momenta 



Fig. 4. 
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ki, ki, respectively. The bare amplitude above is subjected to absorption corrections which 
depend on collision energy and on the spin-parity of the produced central system. 
In the KMR approach the off-diagonal parton distributions are calculated as 

/rv,Q?,^) = R ^ X9[x ^f' y \ , Q , Fit) , (2.11, 

where S1/2 (<2f,/-t 2 ) is a Sudakov-like form factor relevant for the case under consideration 
[24]. The factor R g here is the skewedness parameter and at the RHIC energy the value 
R g ~ 1.4 seems to be relevant. 

In the present calculation we use standard GRV95 collinear gluon distributions [25]. For 
this process we take the renormalization and factorization scales to be ^ 2 R — fi 2 F = M^J 4 as 
in Ref.[22]. Absorption effects are included approximately by multiplying the cross section by 
the gap survival factor Sq = 0.15. More details about EDD production of charm quarks can 
be found in our original paper [22]. In order to compare this predictions with the PHENIX 
data we have applied the same procedure of hadronization and semileptonic decays as in the 
case of inclusive processes described briefly in subsection II A. 

III. NUMERICAL RESULTS 

Let us come now to the presentation of our results. In Fig. 5 we show distribution in 
dielectron invariant mass. We have included several mechanisms. The solid lines represent 
the contribution of charm (upper one) and bottom (lower one) production calculated using 
the Kwiecinski UGDFs and subsequent semileptonic decays which was calculated in the way 
described in subsection II A. The Drell-Yan contribution is shown by the long-dashed line 
and its contribution is comparable to the contribution of semileptonic decays. The gamma- 
gamma contributions (sum of the four contributions of diagrams in Fig. 3) is shown by the 
blue dashed line at the bottom-left corner of the figure. The very small EDD contribution 
is shown for completness by the green dotted line. 

In Fig. 6 and Fig. 7 we present the same distributions but here the calculation of heavy 
quarks is performed with the KMR and Kutak-Stasto UGDFs, respectively. One can see that 
the KMR UGDF gives quite good description of the PHENIX data in the whole considered 
dielectron invariant mass range. Similar results have been obtained with the Kwiecinski 
UGDF except very low dielectron invariant masses. Figure 7 shows that the Kutak-Stasto 
UGDFs reproduce the experimental data only at large dielectron invariant masses but that 
is the region where Drell-Yan mechanism gives significant contribution. As was mentioned 
by the authors of [14], the Kutak-Stasto UGDF is dedicated exclusively to small-x processes 
(x < 10" 2 ), so its use for RHIC is at the border of its applicability, especially for bottom 
quarks. In the calculation of heavy quark/antiquark we have taken m c = 1.5 and = 4.75 
GeV, rather conservative values. 1 

In Fig. 8 we show separately contributions of different photon induced mechanisms shown 
in Fig.3. The amplitude of the double elastic contribution is calculated as explained in 
subsection II C, the other contributions are calculated in the collinear approximation as 
explained in the same subsection. 



Often too small values of heavy quark masses are taken in the calculation to describe the data. 
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M e+e - (GeV) 



FIG. 5: Dielectron invariant mass distribution for proton-proton collisions at yfs = 200 GeV. 
Different contributions are shown separately: semileptonic decay of charm by the upper solid 
line (blue online), semileptonic decay of bottom by the lower solid line (red online), Drell-Yan 
mechanism by the long-dashed line, gamma-gamma processes by the short-dashed line (blue online) 
and the EDD contribution by the dotted line (green online). In this calculation we have included 
azimuthal angle acceptance of the PHENIX detector [1] . 

In Fig. 9 we discuss uncertainties related to the contribution of semileptonic decays. Com- 
plementary the left panel presents uncertainties due to the factorization scale variation as 
described in the figure caption. The right panel shows uncertainties due to the modification 
of the heavy quark masses (m c G (1.25 GeV, 1.75 GeV) and rrib G (4.5 GeV, 5 GeV)). Figure 
10 presents uncertainties of our predictions related to the different models of heavy quark 
fragmentation. Two bands (blue online and red online for charm and bottom, respectively), 
estimated with Peterson fragmentation functions show sensitivity of our results to the vari- 
ation of e parameters in intervals specified in the figure. The long-dashed lines represent 
Braaten et al. perturbative fragmentation model and dotted lines are for the fragmentation 
function proposed by Kartvelishvili et al.. For each of the function we use parameters taken 
from the literature [26, 27]. As one can observe, in comparison with uncertainties discussed 
before there is only a small sensitivity of the results to the fragmentation functions. Some 
small differences start to appear only at large dilepton invariant masses where the error bars 
of the experimental data are realy large. Besides, in the case of bottom quarks such effects 
are almost negligible. 

Since the transverse momenta of electrons can be measured, one can look not only at their 
distributions but also at correlations between them. In Fig. 11 we show two-dimensional dis- 
tribution in transverse momenta of c and c (left panel), D and D mesons (middle panel) and 
e + and e~ (right panel). In contrast to leading-order collinear QCD calculations already 
the distribution at the parton level is dispersed along diagonal. It is further broadened 
by (assumed independent) fragmentation process and even more by (assumed independent) 
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M e+e . (GeV) 

FIG. 6: The same as in Fig. 5 but open charm and bottom components are calculated using KMR 
UGDFs. 




M e+e - (GeV) 



FIG. 7: The same as in Fig. 5 but open charm and bottom components are calculated using Kutak- 
Stasto UGDFs. 

semileptonic decays. So the intial transverse momentum correlations of c and c are prac- 
tically lost when going to the electrons/positrons but are interesting and provide a new 
possibility to test the dynamics of the process and our understanding of QCD at work. 
If the detector can measure both transverse momenta of an electron/positron and its di- 
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M re - (GeV) 



FIG. 8: Contributions of photon-induced mechanisms to the dielectron invariant mass distributions. 
The PHENIX detector limitations are given in the upper-left corner. 




Mrtf < GeV ) M e*e < GeV > 

FIG. 9: The uncertainties of theoretical calculations. The left panel shows the factorization scale 
uncertainties, the lower curve corresponds to ^ 2 F ,^ 2 R = m\ t + m 2 t and the upper curve to fj, R = k 2 , 
Hp = 4mg, where kt is gluon transverse momentum. The right panel shows the quark mass 
uncertainties as indicated in the figure. 

rections, as the STAR detector at RHIC can do, one can construct a distribution in transverse 
momentum of the dielectron pair: pt,sum = Pit + Vit- Our predictions for the semileptonic 
decays and Drell-Yan processes are shown in Fig. 12. Both processes give rather similar 
distributions. To our knowledge the distributions of this type were never measured experi- 
mentally as they cannot easily be compared to the calculations in the collinear approach due 
to its inherent singularities. Obviously this is not the case for the /c t -factorization approach 
discussed in the present analysis. The distribution in pt, S um is not only a consequence of 
gluon transverse momenta, as it is for quark and antiquark production, but invlolves also 
fragmentation process and semileptonic decays. A measurement of this quantity would test 
then all stages of the process. 

With good azimuthal granulation of detectors one could construct distribution in az- 
imuthal angle between electron and positron. Our corresponding predictions are shown in 
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M e+e . (GeV) 

FIG. 10: The uncertainties of theoretical calculations of open charm and bottom related to the 
choice of the fragmentation functions. BCFY means that fragmentation functions from Ref.[17] 
were used. 



charm quarks charm mesons electrons 




P 2t (GeV) p 2 [ (GeV) p y (GeV) 



FIG. 11: Two-dimensional distribution in transverse momenta of cc (left panel), DD (middle panel) 
and e + e~ (right panel). Here Kwiecinski UGDF and Peterson fragmentation function were used. 

Fig. 13. One can see an interesting dependence on the invariant mass of the dielectron pair 
- the smaller the invariant mass the large the decorrelation in azimuthal angle. 

IV. CONCLUSIONS 

In the present analysis we have discussed correlations of charmed mesons and dielectrons 
at the energy of recent RHIC experiments. We have calculated the spectra in dielectron in- 
variant mass, in azimuthal angle between electron and positron as well as for the distribution 
in transverse momentum of the pair. The uncertainties due to the choice of UGDFs, choice 
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FIG. 12: Distribution in transverse momentum of the dielectron pair for semileptonic decays (solid 
line) and Drell-Yan processes (dashed line). Here Kwiecinski UGDF and Peterson fragmentation 
function were used. 



E 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 E 

: p+p at\fs = 200GeV : 




cp e+e . (deg) 



FIG. 13: Distribution in azimuthal angle between electron and positron for semileptonic decays 
(solid line) and Drell-Yan processes (dashed line). Here Kwiecinski UGDFs and Peterson fragmen- 
tation functions were used. 

of the factorization and renormalization scales, choice of the heavy quark masses as well 
as fragmentation functions have been quantified. The uncertainties for UGDFs are larger 
than those for fragmentation functions. We have obtained good description of the dielectron 
invariant mass distribution measured recently by the PHENIX collaboration at RHIC. 

The contribution of electrons from Drell-Yan processes is only slightly smaller than that 
from the semileptonic decays. The distributions in azimuthal angle between electron and 
positron and in the transverse momentum of the dielectron pair from both processes are 
rather similar. We do not find a possibility of a clear separation of both processes. It was 
found that the distribution in azimuthal angle strongly depends on dielectron invariant mass. 
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We have also included exclusive double-diffractive contribution discussed recently in the 
literature. At the rather low RHIC energy it gives, however, a very small contribution to 
the cross section and can be safely ignored. It may not be the case at the LHC energy, as 
the EDD contribution grows much faster than the inclusive cross section. 

The QED double-elastic, double-inelastic, elastic-inelastic and inelastic-elastic processes 
give individually rather small contribution but when added together are not negligible es- 
pecially at low dielectron invariant masses where some strength is clearly missing [1]. 

In the present analysis we have studied correlations between electron and positron and 
in some cases between mesons. It can be also interesting to look at correlations between a 
D meson and electron. This will be a subject of a forthcoming analysis. 



Acknowledgments 

We are indebted to Andre Mischke for exchange of information on recent RHIC results and 
very useful comments. This work was partially supported by the Polish grants MNiSW No. 
N202 2492235 and MNiSW No. N202 236937. 



A. Adare et al. (PHENIX collaboration), Phys. Lett. B670, 313 (2009). 

D. Acosta et al. (CDF II collaboration), Phys. Rev. Lett. 91 (2003) 241804. 
A. Adare et al. (PHENIX collaboration), Phys. Rev. Lett. 97 (2006) 252002. 
M. Cacciari, P. Nason and R. Vogt, Phys. Rev. Lett. 95 (2005) 122001. 

M. Luszczak, R. Maciula and A. Szczurek, Phys. Rev. D79 (2009) 034009. 

M. Luszczak and A. Szczurek, Phys. Rev. D73 (2006) 054028. 

A. Mischke, Phys. Lett. B671, 361 (2009). 

A. Szczurek and G. Slipek, Phys. Rev. D78 (2008) 114007. 

P. Lebiedowicz and A. Szczurek, Phys. Rev. D81 (2010) 036003. 

A. D. Martin, et al., Eur. Phys. J. C39, 155 (2005). 

S. Catani, M. Ciafaloni and F. Hautmann, Nucl. Phys. 366 (1991) 135. 
J. Kwiecihski, Acta Phys. Polon. B33 (2002) 1809. 

M.A. Kimber, A.D. Martin and M. G. Ryskin, Eur. Phys. J. C12, 655 (2000) 

M.A. Kimber, A.D. Martin and M. G. Ryskin, Phys. Rev. D63 (2001) 114027-1. 

K. Kutak and A.M. Stasto Eur. Phys. J. C41, 343 (2005). 

C. Peterson, et. al., Phys. Rev. D27 (1983) 105. 

V.G. Kartvelisbvili, et al., Phys. Lett. B78, 615 (1978). 

E. Braaten, et al., Phys. Rev. D51 (1995) 4819. 

N.E. Adam, et al. (CLEO collaboration), Phys. Rev. Lett. 97 (2006) 251801, 

B. Aubert, et al. (BABAR collaboration), Phys. Rev. D69 (2004) 111104(R). 
Ch-Y. Wong and H. Wang, Phys. Rev. C58 (1998) 376. 

U. D'Alesio and F. Murgia, Phys. Rev. D70 (2004) 074009. 

M. Dress and D. Zeppenfeld, Phys. Rev. D39 (1989) 2536. 

R. Maciula, R. Pasechnik and A. Szczurek, Phys. Lett. B 685, 165 (2010); 

V.A. Khoze, A.D. Martin and M.G. Ryskin, Phys. Lett. B 401, 330 (1997); 



13 



V.A. Khoze, A.D. Martin and M.G. Ryskin, Eur. Phys. J. C 23, 311 (2002). 
[24] A. D. Martin and M. G. Ryskin, Phys. Rev. D 64, 094017 (2001) [arXiv:hep-ph/0107149]. 
[25] M. Gliick, E. Reya and A. Vogt, Z. Phys. C67 (1995) 433. 
[26] W.M. Yao, et al. [Particle Data Group], J. Phys. G33 (2006) 1. 
[27] J. Baines, et al., HERA-LHC Workshop Proceedings, 

Heavy Quarks: Summary Report, arXiv:hep-ph/0601164v2. 



14 



